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ABSTRACT 


Thia  report  summarizes  the  results  of  a study  on  various  nonlinear 
phenomena  In  high-temperature  plasmas;  in  particular,  explosive  instabil- 
ities wave— energy  transfer  in  nonlinear  plasmas,  turbulence,  numerical 

solution  of  nonlinear  equations  in  plasma  physics,  plasma  heating  by  neutral 
injection,  feedback  stabilization  and  control  of  high-temperature  plasmas. 


1.  INTRODUCTION 

Recent  experiments  in  plasma  physics  involving  the  Interaction  of 
high-intensity  electromagnetic  radiation  with  plasmas,  heating  of  plasmas 
by  Injecting  high-energy  particles,  and  high-current  toroidal  plasmas 
often  produced  results  which  cannot  be  explained  using  linear  theory.  A 
typical  example  is  the  explosive  instabilities  resulting  from  the  interac- 
tion of  a system  of  positive  and  negative  energy  waves  in  a plasma.  The 
Induced  emission  of  hlgh-lntenslty  electromagnetic  radiation  from  a high- 
temperature  plasma  is  not  predictable  by  means  of  linear  models.  Also, 
the  onset  of  turbulence  in  a plasma  is  basically  a nonlinear  phenomena. 

It  is  evident  that  nonlinear  phenomena  will  play  an  increasingly  Important 
role  in  the  future  development  of  both  theoretical  and  applied  plasma  phy- 
sics. 

From  January  1,1974  through  December  31,1978,  studies  of  nonlinear 
phenomena  in  high-temperature  plasmas  under  the  support  of  AFOSR  grant  No. 
74-2662  have  been  concentrated  in  the  following  three  main  areas:  (a)  non- 
linear plasma  theory,  (b)  plasma  heating  by  neutral-particle  Injection  and 
(c)  feedback  stabilization  and  control  of  plasmas.  The  results  in  each 
of  these  areas  will  be  sumswrlzed  categorically  below. 

2.  NONLINEAR  PLASMA  THEORY 

Studies  in  nonlinear  plasmas  began  with  an  investigation  of  nonlinear 
wave-wave  Interactions  in  plasmas  with  particular  emphasis  on  explosive 
instabilities.  This  work  was  motivated  from  a number  of  experiments  on 
the  wave-heating  of  plasmas  in  which  wave  amplitudes  with  nonexponential 
growth  rates  were  observed.  Later,  this  work  was  broadened  to  Include 
plasma  turbulence  theory  from  the  viewpoint  of  bifurcation  theory  and  soli- 
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ton  Interactions.  In  parallel  with  this  work,  attempts  were  made  In 
developing  efficient  accurate  numerical  methods  for  solving  the  nonlinear 


plasma  equationa. 

2.1  Explosive  Instabilities  and  Wave-Energy  Transfer  In  Nonlinear  Plasmas 
It  Is  well-known  that  In  nondlssipative  plasmas  with  nonlinear ly  Inter- 
acting positive  and  negative  energy  waves,  the  wave  amplitudes  could  grow 
at  rates  faster  than  exponential  [l]-[3].  Moreover,  some  of  the  waves  may 
exhibit  an  explosive  behavior  In  the  sense  that  their  wave  amplitudes  tend 
to  Infinity  In  a finite  amount  of  time.  Of  course.  In  physical  situations, 
various  saturation  phenomena  take  place  when  the  wave  amplitude  becomes  suf- 
ficiently large.  In  plasmas,  explosive  instabilities  have  been  observed 
In  a number  of  experiments  [4], [5].  During  recent  years,  theoretical  studies 
In  this  area  have  been  focused  on  the  derivation  of  necessary  and/or  suf- 
ficient conditions  for  the  existence  or  nonexistence  of  explosive  Instabili- 
ties In  plasmas  with  nonlinear ly  interacting  waves.  For  the  case  of  non- 
diaslpatlve  nonlinear  plasmas  having  waves  with  deterministic  phases,  con- 
ditions for  the  existence  of  explosive  Instabilities  have  been  established 
for  the  case  with  weakly  nonlinear  wave-wave  Interactions  [6].  Wllhelmsson 
at  al  [7]  derived  a necessary  condition  for  explosive  instability  for  three- 
wave  interactiona  taking  Into  account  the  effect  of  linear  damping  or  growth. 
In  a study  supported  by  AFOSR  grant  72-2303,  Wang  [8]  obtained  a sufficient 
condition  for  the  nonexistence  of  explosive  Instabilities  by  first  deriving 
appropriate  bounds  for  the  solutions  of  the  differential  equations  govern- 
ing the  dynamics  of  nonlinear  wave-wave  Interactions.  In  spite  of  the 
large  amount  of  existing  work  pertaining  to  the  explosive  Instabilities  In 
various  types  of  nonlinear  plasmas,  no  attempt  has  been  made  to  unify  the 
existing  results  or  to  approach  the  problem  from  a general  viewpoint.  Most 


works  dealt  with  very  special  models  of  nonlinear  plasmas  and  classical 
mathematical  techniques  were  used  In  the  analysis. 

In  this  study,  we  tried  to  approach  the  problem  from  the  standpoint 
of  qualitative  theory  of  differential  equations.  The  objective  Is  to  ob- 
tain some  general  results  pertaining  to  the  existence  and  nonexistence  of 
explosive  solutions  of  those  classes  of  differential  equations  arising  In 
nonlinear  plasma  theory.  Then,  we  proceeded  to  Investigate  perticular 
mathematical  models  of  nonlinear  plasmas  in  the  light  of  the  general  re- 
sults. 

The  starting  point  of  this  study  Is  based  on  the  observation  that  the 
differential  equations  describing  almost  all  the  nonlinear  plasmas  can  be 
formulated  In  the  form  of  a nonlinear  evolution  equation  of  the  form: 

g « Ax  + f (x),  (1) 

defined  on  some  suitable  Hilbert  space  H,  where  A Is  a linear  operator  on 

H into  H,  and  f Is  a nonlinear  operator  having  the  homogeneity  property 
v 

f(cx)  • c f(x)  for  some  k > 0.  For  example.  In  the  case  of  a system  of 
a finite  number  of  nonlinear ly  Interacting  waves  In  a plasma,  equation  (1) 
corresponds  to  a system  of  ordinary  differential  equations  defined  on  the 
n-dlmens lonal  complex  vector  space  Cn  with  x being  the  complex  wave-ampli- 
tude vector.  On  the  other  hand,  for  a nonlinear  plasma  described  by  the 
nonlinear  fluid  equations,  (1)  corresponds  to  a system  of  partial  differen- 
tial equations  In  which  A is  a linear  spatial  differential  operator  and  f 
la  a nonlinear  homogeneous  spatial  differential  operator  with  degree  of 
homogeneity  k - 2. 

In  this  study,  attention  was  focused  on  developing  general  sufficient 
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conditions  for  the  existence  and  nonexistence  of  explosive  solutions  of 
(1).  Results  were  obtained  by  projecting  the  solutions  of  (1)  onto  the 
unit  sphere  in  the  space  H and  then  studying  the  behavior  of  the  solutions 
which  lie  along  rays  emanating  from  the  origin  of  the  space.  It  was  shown 
that  explicit  conditions  for  which  the  ray  solutions  are  explosive  or  non- 
explosive  can  be  obtained  by  solving  an  eigenvalue  problem  associated  with 
the  operators  A and  f.  When  A and  f belong  to  certain  special  classes  of 
operators,  in  particular,  potential  operators  and  compact  operators,  the 
existence  of  explosive  solutions  can  be  established  directly  in  terms  of 
the  parameters  of  (1).  Moreover,  the  explosion  times  for  the  solutions 
^ can  be  computed  exactly.  The  aforementioned  classes  of  operators  arise  in 

nonlinear  wave-wave  interactions  in  plasmas.  The  details  of  these  results 
are  described  in  paper  (P-2). 

In  order  to  determine  the  implications  of  the  general  theoretical  re- 
sults, studies  were  made  on  the  nonlinear  interaction  of  a finite  number 
of  waves  in  plasmas.  It  was  shown  that  the  wave  energy  transfer  along 
any  ray  solution  is  always  unidirectional  or  irreversible,  l.e.  a wave 
gains  or  loses  energy  along  a ray  solution  at  all  times.  This  property  is 
significant  from  the  application  standpoint,  since  it  Implies  that  the  ray 
solutions  represent  efficient  paths  for  transferring  energy  from  one  wave 
to  another.  Evidently,  this  type  of  solutions  may  be  selected  for  heat- 
ing a plasma  to  high  temperatures  (for  example,  by  transferring  energy 
from  an  Intense  external  wave  to  the  plasma  waves).  Detailed  results 
were  worked  out  for  the  special  case  of  nonlinear  two-wave  and  three-wave 
Interactions.  The  significance  of  these  results  for  the  nonlinear  re- 
sonant interaction  of  three  longitudinal  electromagnetic  waves  in  a magne- 
tised collisionless  plasma  was  also  studied.  The  results  of  this  aspect 


of  study  were  published  in  paper  (P-5).  Pursuing  further  along  these  di- 
rections, it  was  found  during  the  first  part  of  1975  that  the  unidirectional 
energy  transfer  property  not  only  holds  along  the  ray  solutions  of  a general 
wave-wave  interacting  system,  but  also  along  those  solutions  whose  phases 
are  invariant  with  time.  This  suggests  that  in  the  wave-heating  of  plasmas, 
the  initial  phases  of  the  waves  could  be  chosen  in  such  a way  to  achieve 
more  efficient  transfer  of  energy.  The  results  of  this  study  were  first 

a 

presented  at  the  Second  International  Congress  on  Waves  and  Instabilities 
in  Plasmas  held  in  Innsbruck  (Paper  (CP-1))  and  published  later  in  Physica 
(Paper  (P-7)).  Also,  the  results  pertaining  to  the  wave  energy  transfer 
in  time-varying  nonlinear  wave-wave  interacting  systems  were  obtained  and 
published  in  paper  (P-1). 

In  late  1975,  an  attempt  was  made  to  derive  more  general  sufficient 
conditions  for  the  existence  or  nonexistence  of  explosive  instabilities. 

In  early  1976,  we  obtained  a set  of  general  conditions  for  stability  and 
for  explosive  instability  by  using  certain  relations  between  an  abstract 
function  and  its  Gateaux  differential.  It  was  shown  that  these  results 
are  applicable  to  a large  class  of  equations  arising  from  physical  systems 
involving  nonlinear  wave  propagation,  diffusion  or  transport  delays.  A 
special  case  of  these  results  is  the  well-known  Krasovskll's  stability 
theorem  for  finite-dimensional  ordinary  differential  equations.  These  re- 
sults were  published  in  the  paper  (P-8). 

In  the  study  of  explosive  Instabilities  and  wave-energy  transfer  in 
specific  wave-wave  interacting  systems,  results  were  obtained  for  the  para- 
metric decay  of  an  intense,  coherent  electromagnetic  wave  into  a plasma 
wave  and  scattered  electromagnetic  waves  in  a homogeneous  plasma  (stimulated 
Raman  scattering).  This  study  includes  the  coupling  of  secondary  waves 
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consisting  of  s plasma  wave  harmonic  and  Its  associated  mixed  electromagne- 
tic-electrostatic sidebands.  Explicit  expressions  were  obtained  for  the 
wave  amplitudes  and  nonlinear  growth  rates  when  t > 1/y,  where  Y Is  the 
linear  growth  rate  of  the  fundamental.  It  was  shown  that.  In  this  regime, 
the  nonlinear  growth  rate  Is  2y/3.  The  results  are  described  In  report 
CR-2). 

Finally,  a general  formalism  was  developed  to  describe  the  nonlinear 
evolution  associated  with  the  parametric  decay  of  an  Intense,  coherent  elec- 
tromagnetic wave  Into  an  electrostatic  wave.  Its  second  harmonic,  and  scat- 
tered electromagnetic  waves  in  a homogeneous  plasma.  Neglecting  the  ef- 
fects of  pump  depletion  and  assuming  all  waves  are  coherent,  two  classes  of 
solutions  were  found.  One  class  consists  of  explosively  unstable  solutions 
and  the  other  consists  of  growing  aperiodic  oscillatory  solutions.  The 
details  are  given  In  report  (R-3) . These  results  have  applications  In 
ionospheric  modification  experiments  and  In  the  interaction  of  laser  beams 
with  plasmas. 

2.2  Plasma  Turbulence 

A possible  approsch  to  developing  a theory  for  strongly  turbulent 
plasmas  is  via  the  study  of  the  Interaction  of  solltons.  In  1973,  Kaplan 
and  Tsytovlch  [7]  showed  by  means  of  quasl-llnear  theory  that  the  energy 
in  a Langmuir-turbulent  plasms  tends  to  accumulate  at  very  small  wave  num- 
bers. However,  In  view  of  the  results  of  Vedenov  and  Rudako  [8],  It  Is 
apparent  that  this  nearly  uniform  plasma  state  Is  unstable.  In  fact,  mod- 
ulatlonal  Instability  can  occur  which  leads  to  local  density  depressions 
along  with  Increases  In  the  energy  density  of  the  Langmuir  waves.  In  1972, 
Zakharov  [9]  formulated  the  basic  equations  governing  the  aodulatlonal 
instability  in  which  the  high-frequency  Langmuir  oscillations  are  coupled 
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to  the  Ion  motion  by  the  ponderomotive  force.  It  was  found,  in  the  one- 
dimensional  case,  that  these  equations  possess  soliton  solutions  in  which 
the  Langmuir  oscillations  are  trapped  in  regions  of  lower  density.  Later, 
Kingsep  et  al  [10]  suggested  that  strong  Langmuir  turbulence  could  be  des- 
cribed in  terms  of  a system  of  interacting  Langmuir  solitons.  A consider- 
able amount  of  work  both  analytical  and  numerical,  has  been  done  based  on 
this  idea  for  the  one-dimensional  case  [11]-[13].  It  was  hoped  that  simi- 
lar theory  could  be  developed  for  the  three-dimensional  case.  During  1977, 


we  obtained  general  sufficient  conditions  for  the  existence  or  nonexistence 


of  a class  of  multi-dimensional  Langmuir  solitons  based  on  a generalized 
version  of  the  Zakharov's  model.  It  was  shown  that  these  solitons  having 
properties  similar  to  those  of  the  one-dimensional  case  could  exist  under 
subsonic  conditions.  Also,  periodic  nonlinear  travelling  waves  could  form 
under  both  subsonic  and  supersonic  conditions.  The  results  were  applied 
to  the  special  cases  with  the  usual  ponderomotive  force  and  with  ion  density 
saturation.  The  results  were  published  in  paper  (P-9). 

In  early  1978,  we  began  exploring  a new  direction  for  developing  a de- 
terministic theory  for  plasma  turbulence.  This  work  was  motivated  from 
Ruelle-Takens*  proposal  [14]  on  the  mathematical  characterization  of  fluid 
turbulence  in  terms  of  "strange-attractor  solutions".  This  class  of  solu- 
tions corresponds  to  those  trajectories  in  the  state  space  of  the  dynamical 
model  which  are  attracted  to  a nonempty  set  ("strange  attractor").  This 
set  is  neither  an  equilibrium  set  nor  a periodic  orbit.  Moreover,  the  auto- 
correlation function  p(T)  for  these  solutions  tends  to  zero  as  the  separa- 
tion time  T -*■  «.  In  fluid  turbulence,  Lorenz's  simplified  model  for  ther- 
mal convection  In  the  atmosphere  seems  to  have  solutions  exhibiting  such 
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elasticity.  It  has  the  capability  of  achieving  high  accuracy  with  a small 
number  of  approximating  equations.  Here,  the  numerical  solution  of  the 
Vlasov-Poisson  system  of  nonlinear  equations  for  a multispecies  plasma  by 
means  of  the  finite-element  Galerkin  approximation  was  considered.  The  fol- 
lowing numerical  algorithm  was  proposed:  First,  the  time  Interval  Is  parti- 
tioned Into  subintervals.  At  the  beginning  of  each  subinterval,  the  Poisson's 
equation  with  the  charge  density  at  the  end  of  the  previous  subinterval  Is 
solved.  Then,  the  linearized  Vlasov's  equations  are  solved  using  the  finite- 
element  Galerkin  approximation.  Explicit  error  estimates  were  derived. 

Using  these  estimates.  It  Is  shown  that  the  approximate  solutions  generated 
by  the  proposed  algorithm  converge  to  the  true  solution.  The  rate  of  conver- 
gence was  also  established.  To  gain  some  computational  experience  with  the 
proposed  algorithm,  a computer  code  was  developed  for  the  case  of  an  electron- 
ion  collisionless  plasma.  The  numerical  results  compared  favorably  with  the 
theoretical  error  estimates.  The  details  are  presented  in  the  Ph.D.  thesis 
by  Ezzuddln  (T-l) . 

A study  was  also  made  on  the  numerical  solution  of  the  equations  for 
non-equilibrium  nozzle  flow  In  CO2-N2  gasdynamlc  lasers.  A computational 
method  based  on  optimal  control  theory  was  developed.  In  this  method,  the 
nozzle  area  and  the  gradient  of  the  gas  temperature  are  treated  as  the  state 
and  control  variables  respectively,  and  the  non-equillbrlum  effects  due  to 
the  kinetic  processes  are  considered  In  both  the  upstream  and  downstream 
of  the  nozzle  throat.  This  approach  bypasses  the  difficulties  due  to  the 
sonic-point  singularity.  The  results  of  this  study  are  given  in  the  M.S. 
thesis  by  Kasul  (T-2). 
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3.  PLASMA  HEATING  BY  NEUTRAL  INJECTION 

During  Che  first  hslf  of  1975,  the  principal  Investigator  visited 
the  Department  of  Plasma  Physics  and  Controlled  Fusion,  Center  of  Nuclear 
Studies  at  Fontenay-aux-Rosea , France.  A study  In  collaboration  with 
M.  Cotsaftls  was  made  on  the  optimal  heating  of  plasmas  by  neutral  injec- 
tion. A mathematical  model  for  the  time  evolution  of  the  spatially  averaged 
electron  and  ion  temperatures  was  used  in  this  study.  The  model  takes  Into 
account  the  major  sources  and  sinks  which  are  Important  in  the  KeV  range. 
First,  the  dynamics  of  plasma  heating  was  analyzed  in  the  electron-ion 
temperature  plane  for  the  case  with  constant  energy  and  current  neutral 
injection.  By  investigating  the  behavior  of  the  right-hand-sides  of  the 
temperature  equations,  the  exact  regions  of  heating  and  cooling  of  electrons 
or  ions  were  found.  It  was  shown  that  there  exists  a region  In  which  heat- 
ing of  both  the  electrons  and  ions  takes  place,  and  the  size  of  this  region 
Increases  with  the  injection  current  or  energy.  Figure  1 shows  the  qualita- 
tive behavior  of  typical  trajectories  in  the  electron- ion  temperature  plane 
for  the  case  with  constant  current  and  energy  injection. 

The  second  phase  of  this  study  pertains  to  various  problems  associated 
with  the  optimization  of  the  neutral  injection  parameters.  In  particular, 
the  problem  of  minimizing  the  total  input  energy  to  achieve  the  highest  ion 
temperatures  at  the  end  of  a specified  injection  period  was  studied.  The 
problem  of  time-optimal  heating  was  also  Investigated.  Here,  an  Injection 
current  and/or  energy  program  was  sought  to  attain  a given  value  of  ion- 
tesperature  in  a »<«■*■»«»■  amount  of  time  subject  to  certain  limitations  on 
the  —rtimmi  values  of  the  injection  current  and  energy.  These  problems  were 
studied  in  the  spirit  of  optimal  control  theory.  The  results  of  this  study 
were  applied  to  the  Tokamak  device  at  Fontenay-aux-Roses.  Good  agreement 
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was  obtained  between  the  theoretical  and  experimental  results.  These  results 
were  presented  at  the  American  Physical  Society,  Plasma  Physics  Division 
Meeting  In  St.  Petersburg,  Florida,  In  1975  (Conference  Papers  CPA-3  and 
CPA-4).  They  are  also  described  In  the  report  (R-l). 

This  study  was  continued  during  M.  Cotsaftls'  visit  to  U.C.L.A.  In 
1976.  The  main  Idea  of  the  continuing  study  was  motivated  from  the  results 
of  the  neutral  Injection  experiments  at  Fontenay-aux-Roses  during  1976.  It 
was  observed  that  the  electron  temperature  did  not  vary  appreciably  during 
the  neutral  Injection  period  while  the  Ion  temperature  Increased  substan- 
tially. This  could  be  attributed  to  the  fact  that,  when  the  plasma  current 
was  sufficiently  large,  the  electrons  were  In  a high- loss  regime  dominated 
by  trapped-electron  Instability.  In  this  esse,  the  two-temperature  model 
developed  earlier  could  be  simplified  by  setting  the  electron  temperature  at 
Its  equilibrium  value.  Based  on  this  simplified  model,  we  solved  completely 
the  problem  of  finding  an  Injection  current  program  1*1(0  satisfying  an 
amplitude  constraint  such  that  the  Ions  are  heated  to  a specified  temperature 
while  minimizing  the  heating  time  or  the  Injection  energy  over  a given  time 
Interval.  The  dynamics  and  stability  of  an  lon-temperature  feedback  control 
system  were  also  analyzed  for  various  forms  of  nonlinear  feedback  controls. 

The  analysis  Included  the  effect  of  measurement  time-delay  on  system  stability. 
The  results  suggest  the  possibility  of  using  automatic  control  for  close  regu- 
lation of  the  ion- temperature  by  means  of  neutral  Injection.  The  problem  of 
Identifying  the  Ion  diffusion  law  by  measuring  the  Ion  temperature  over  a 
finite  period  of  time  was  also  explored.  The  results  of  this  study  are  de- 
scribed in  paper  (P-10)  and  were  presented  at  the  1976  Annual  Meeting  of  the 
Plasma  Physics  Division  of  the  American  Physical  Society  In  San  Francisco. 
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Finally,  tha  problem  of  optimal  timing  and  currant  programming  of 
neutral  beam  Injection  to  achieve  the  most  efficient  heating  of  plaamaa 
vaa  investigated.  It  was  shown  that  this  problem  can  be  formulated  as  a 
multi-stage  optimal  control  problem.  First,  certain  general  results  per- 
taining to  the  existence  and  characterization  of  optimal  controls  were 
obtained.  Then,  they  were  applied  to  the  optimal  neutral  Injection  heating 
problem.  The  results  are  presented  In  the  Ph.D.  dissertation  by  K.  Tomi- 
yama  (T-3). 

4.  FEEDBACK  STABILIZATION  AND  CONTROL  OF  HIGH-TEMPERATURE  PLASMAS 

In  the  production  of  hlgh-temperature  laboratory  plasmas.  It  Is 
desirable  to  be  able  to  vary  or  regulate  various  plasma  parameters  such  as 
temperature,  density,  etc.,  as  required  by  a particular  experiment.  Due  to 
the  Inherent  Instabilities  and  physical  constraints,  the  ranges  of  allowable 
parameter  variations  are  narrow.  Consequently,  such  plasmas  cannot  be  used 
for  a wide  range  of  experiments.  By  Introducing  feedback  controls  to  regu- 
late the  plasma  parameters.  It  may  be  possible  to  extend  the  useful  parameter 
range,  thus  enhancing  the  usefulness  of  the  laboratory  plasma. 

The  main  objective  of  this  study  Is  to  develop  a theory  for  the  feed- 
back stabilization  and  control  of  hlgh-temperature  plasmas  from  a unified 
control  theoretic  viewpoint.  The  specific  areas  of  this  study  will  be  dis- 
cussed separately  below. 

4.1  Feedback  Stabilization 

Given  a mathematical  model  for  a plasma  to  be  stabilized  by  means  of 
some  form  of  external  controls  (for  example,  maripulatable  currents,  electric 
fields,  etc.),  the  following  basic  questions  should  be  answered  before  attempt- 
lng  to  Implement  any  feedback  controls:  (1)  Do  the  controls  affect  all  the 

13 


■ 


node.  It  is  desirable  to  have  a control  system  such  that  each  control 
(input  conaand)  affects  only  one  mode  and  each  mode  is  affected  by  only 
one  control.  Thus,  there  is  a one-to-one  correspondence  between  the  input 
commend  and  the  modes  to  be  controlled.  This  may  be  accomplished  by  intro- 
ducing suitable  feedback  signals  within  the  plasma  system.  During  the 
period  covered  by  this  report,  a study  was  made  on  the  noninteracting  con- 
trols for  high-temperature  plasmas.  Attention  was  focused  on  developing 
theories  for  noninteracting  modal  controls  for  plasmas  describable  by  the 
MHD  equations  and  the  Vlasov  equation.  Both  the  MHD  modes  and  drift  modes 
ware  considered.  Results  were  obtained  pertaining  to  the  following  basic 
questions:  (1)  Given  a mathematical  model  for  the  plasma,  under  what  condi- 
tions can  one  achieve  the  noninteracting  controls?  (2)  What  should  be  the 
form  end  the  location  of  the  controls  and  sensors  for  implementing  the  non- 
interacting controls?  The  details  are  given  In  (C-l)  and  were  presented  at 
the  First  IEEE  International  Conference  on  Plasma  Science  (see  Conference 
paper  CPA-1). 

The  practical  significance  of  these  results  is  that,  once  the  non- 
Intaractlng  modal  controls  have  been  achieved,  the  problem  of  stabilisation 
of  unstable  modes  becomes  straightforward.  Moreover,  each  mode  can  be  sta- 
bilised or  enhanced  separately  by  introducing  feedback  of  the  modal  ampli- 
tude of  that  mode  only,  without  affecting  the  remaining  modes.  This  festure 
is  very  useful  from  an  experimentalist's  viewpoint,  since  each  mode  can  be 
completely  controlled  and  studied  in  detail  at  the  experimentalist's  dis- 
cretion. In  plasma  experlmsnts,  the  noninteraction  controls  are  also  useful 
In  achieving  independent  adjustment  of  the  Important  plasma  parameters. 
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4.2  Optimal  Control 


During  th«  period  covered  by  this  report,  various  new  optimal  control 
problems  motivated  by  plasma  control  were  studied.  In  particular,  the  prob- 
lem of  optimal  control  of  systems  governed  by  parabolic  differential  equa- 
tions with  boundary  conditions  Involving  time-delays  was  Investigated.  This 
problem  was  derived  from  that  of  confining  a collision-dominated  plasma  by 
means  of  reflective  barriers  (such  as  a magnetic  mirror  field)  at  the  plasma- 
domain  boundary.  The  details  are  published  In  paper  (P-3).  In  the  same  vein 
optimal  control  problems  for  time-lag  systems  with  time-lag  controls  were 
studied  (see  papers  (P-4)  and  (P6)).  Finally,  it  was  shown  that  the  problem 
of  optimal  confinement  of  plasmas  by  means  of  external  electromagnetic  fields 
can  be  formulated  mathematically  as  a problem  in  the  optimization  of  set 
functions;  i.e.,  by  choosing  a geometric  domain  In  a certain  class  to  mini- 
mize a given  functional.  General  necessary  and  sufficient  conditions  for 
optimality  were  obtained  for  this  class  of  problems  (see  (P-11),  (CP-3)  and 
(T-4)  for  details). 
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4.3  Parameter  Identification  In  High-Temperature  Plasmas 

In  the  generation  of  high-temperature  plasmas.  It  Is  of  Importance  to 
estimate  the  basic  plasma  paramsters,  such  as  density  and  temperature,  from 
the  physically  observable  quantities.  In  the  implementation  of  feedback  con- 
trols for  plasmas,  it  Is  also  necessary  to  have  good  estimates  of  the  Impor- 
tant plasma  paramsters,  which  are  usually  la  the  nonequll ibr lum  state. 

There  are  two  basic  aspects  to  this  problem,  namely,  devising  suitable  meas- 


urement or  diagnostic  techniques  for  plasmas,  and  processing  the  experimental 
data  to  obtain  the  best  estimates  of  the  plasma  parameters.  Existing  work 
In  plasma  diagnostics  deals  primarily  with  developing  new  measurement  tech- 
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